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We investigate the origin of diffusion in non-chaotic systems. As an example, we consider 1-d map 
models whose slope is everywhere 1 (therefore the Lyapunov exponent is zero) but with random 
quenched discontinuities and quasi-periodic forcing. The models are constructed as non-chaotic 
approximations of chaotic maps showing deterministic diffusion, and represent one-dimensional ver- 
sions of a Lorentz gas with polygonal obstacles (e.g., the Ehrenfest wind tree model). In particular, 
a simple construction shows that these maps define non-chaotic billiards in space-time. The models 
exhibit, in a wide range of the parameters, the same diffusive behavior of the corresponding chaotic 
versions. We present evidence of two sufficient ingredients for diffusive behavior in one-dimensional, 
non-chaotic systems: i) a finite-size, algebraic instability mechanism, and ii) a mechanism that 
suppresses periodic orbits. 

PACS numbers: 05.45.-a, 05.60.-k 
Keywords: Diffusion, Chaos 



I. INTRODUCTION 



The presence of randomness in many natural processes 
is traditionally ascribed to the lack of perfect control on 
observations and experiments. Then, usually, "disorder" 
and "irregular behaviors" are considered only as unavoid- 
able factors which, in phenomena modellizations, can 
be described effectively through external sources (ran- 
dom numbers in computer simulations, noise in stochas- 
tic equations etc.). However randomness may also have 
another origin associated to local instabilities of the de- 
terministic microscopic dynamics of a system, i.e. micro- 
scopic chaos (MC). For instance, macroscopic diffusion of 
a particle in a fluid can be interpreted as a consequence 
of chaos on microscopic scales generated by collisions of 
the particle with the surrounding others. 

The possibility of describing the Brownian behavior 
in terms of deterministic motion has stimulated a large 
amount of work in the field of deterministic diffusion 
to study purely deterministic systems exhibiting asymp- 
totically a linear growth of the mean square displace- 
ment UJ- L3- La] • The deterministic interpretation of trans- 
port phenomena establishes a close connection between 
transport and chaos theory, and, at least in principle, 
macroscopic properties such as transport coefficients (vis- 
cosity, thermal and electrical conductivity, diffusion etc.) 
could be directly computed through indicators of chaos 
[3- EL El or by applying the standard techniques of chaos 
theory such as periodic orbit expansion |7ll8ll9lll0l|. How- 
ever, considering randomness and transport in general as 
signatures of MC raises subtle and fundamental issues in 



statistical physics that have to be carefully discussed an 
investigated. First, there is the main problem of conceiv- 
ing and realizing experiments able to prove, beyond any 
doubt, the existence of MC in real systems. An experi- 
ment in this direction was attempted in ref. but the 
conclusions were object of criticism by many authors, 
as we discuss later. Besides, for real systems, the def- 
inition of MC appears to be ambiguous. Indeed, real 
systems involve a large number of degrees of freedom 
and their model description entails taking the thermody- 
namic limit but, in this limit, the concept of Lyapunov 
exponent becomes metric dependent so not well estab- 
lished Moreover, the large number of degrees of free- 
dom poses severe limitations also in the applicability of 
standard techniques to dis ting uish between chaotic and 
stochastic signals 0, 0, when data are extracted 
from real experiments. In fact, the noise-chaos distinc- 
tion requires a extremely fine resolutions in the obser- 
vations that is impossible to reach when one deals with 
systems with large dimensionality 0, 0, 0] . 
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References ^2 113 summarized an ingenious experi- 
ment devised to infer the presence of chaos on micro- 
scopic scales. In the experiment the position of a Brown- 
ian particle in a fluid was recorded at regular time inter- 
vals. The time series were processed by using the same 
techniques developed for detecting chaos from data anal- 
ysis 0, . Authors measured a positive lower bound 
for the Kolmogorov-Sinai entropy, claiming that this was 
an experimental evidence for the presence of microscopic 
chaos. However in the papers 0] and it was ar- 
gued that a similar result holds for a non-chaotic system 
too, hence the experiment cannot provide a conclusive 
evidence neither for the existence of microscopic chaos 
nor for the relevance of chaotic behavior in diffusion phe- 
nomena. Such papers shown that, through the analysis 
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performed in Refs. [Til Il9j . there is no chance to ob- 
serve differences in the diffusive behavior between a gen- 
uine deterministic chaotic systems as the 2-d Lorentz gas 
with circular obstacles [(| and its non chaotic variant: the 
wind-tree Ehrenfest model. The Ehrenfest model consists 
of free moving independent particles (wind) that scatter 
against square obstacles (trees) randomly distributed in 
the plane but with fixed orientation. Due to collisions, 
particles undergo diffusion, however their motion can- 
not be chaotic because a reflection by the flat planes of 
obstacles does not produce exponential trajectory sep- 
aration; the divergence is at most algebraic leading to 
zero Lyapunov exponent. Such considerations can be 
extended to squares obstacles with random orientations 
and to every polygonal scatterers, so there is a whole 
class of models where diffusion occurs in the absence of 
chaotic motion. A difference between the 2-d chaotic and 
non-chaotic Lorentz gas relies on the presence of periodic 
orbits [171 ]. 

The question that naturally arises "What is the mi- 
croscopic origin of the diffusive behavior?" remains still 
open. In this paper, focusing on this problem, we extend 
the analysis of Dettman and Cohen showing that stan- 
dard diffusion may also occur in one-dimensional models 
where every chaotic effect is surely absent. These sys- 
tems, being one-dimensional, can be much easily stud- 
ied than wind-tree models without spoiling the essence 
of the problem. Accurate results are obtained with rela- 
tively small computational efforts, and the interpretation 
of the results is clearer since we are able to better control 
and quantify the effects of the local instabilities. We shall 
see that, in such models, the relevant ingredients to ob- 
tain diffusion in the absence of deterministic chaos is the 
combined action of a "finite size instability", quenched 
disorder and a quasi-periodic perturbation. Where by 
"finite size instability" we mean that infinitesimal pertur- 
bations are stable, while perturbations of finite size can 
grow algebraically. The interplay between these three 
features guarantees the system to be still diffusive as it 
would be chaotic. This outcome is another striking effect 
of finite size instabilities which make the behaviors of cer- 
tain non-chaotic system similar, to some extent, to those 
generated by genuine chaotic systems ( H3 El IH H| ) . 
In ref. [24[ the reader can find a study of an interesting 
system with weak mixing properties and zero Lyapunov 
exponent, i.e. an intermediate case between integrability 
and strong mixing. 

The paper is organized as follows. In section II, we first 
recall the general features of 1-d models showing deter- 
ministic diffusion. Then we present two non-chaotic, one- 
dimensional models that exhibit deterministic diffusion 
when properly perturbed. In section III we characterize 
numerically the diffusion properties of the models and 
their dependence on the external perturbation. Section 
IV shows how these models can be viewed as non-chaotic 
space-time billiards, and discusses in more detail the role 
of periodic orbits. Section V is devoted to conclusions 
and remarks. 



II. NON-CHAOTIC MODELS FOR 
DETERMINISTIC DIFFUSION 

Motivated by the problem of deterministic diffusion in 
non-chaotic Lorentz systems (i.e. with polygonal obsta- 
cles) we introduce models that, somehow, represent their 
one-dimensional analog. We begin considering one of the 
simplest chaotic model that generates deterministic dif- 
fusion: a 1-d discrete-time dynamical system on the real 
axis 



+ F(x t -[x t ]), 



(1) 



where x* is the variable performing the diffusion (position 
of a point-like particle) and [. . .] denotes the integer part. 
F(u) is a map defined on the interval [0, 1] that fulfills 
the following properties 

i) |F'(u)| > 1, so the system has a positive Lyapunov 

exponent. 

ii) F(u) must be larger than 1 and smaller than for 

some values of it, so there exists a non vanishing 
probability to escape from each unit cell (a unit 
cell of real axis is every interval Ct = [£, I + 1] , with 

tez). 

iii) F r (u) = 1 — F/(l — u), where Fi and F r define the map 

in u E [0, l/2[ and u £ [1/2,1] respectively. This 
anti-symmetry condition with respect to u = 1/2 is 
introduced to prevent the presence of a net drift. 

The map, u t+1 = F(u l ) (mod 1) is assumed to be also 
ergodic. One simple choice of F is 



F(u) 



2(1 + a)u if ue [0,1/2] 

2(1 + a)(u - 1) + 1 if uE [1/2,1] 



(2) 



where a > is the parameter controlling the instability. 
This model has been so widely studied, both analytically 
and numerically 3, 25, 26( that can be considered as a ref- 
erence model for deterministic diffusion. The determin- 
istic chaos which governs the evolution inside each unit 
cell is responsible of the irregular jumping between cells 
and the overall dynamics appears diffusive, in the sense 
that asymptotically the variance of x scales linearly with 
time. Following the literature on deterministic diffusion 
we consider systems with no net drift (i.e. (x )/t — > for 
t — > oo). In this paper we will see that there is no need 
to invoke chaos to induce diffusive behavior in 1-d maps, 
as already shown in ^3] m the context of 2-d billiards. 

The construction of the models begins by dividing the 
first half of each unit cell Ci into N micro-intervals [£ + 
S,n-lJ + Cn[ , n = 1,...,N, with Co = < Ci < 6 < 
• • • < £jv— i < £jv = 1/2- ln each micro-interval the maps 
are defined by prescribing a linear function with unit 
slope. The map in the second half of the unit cell is 
determined by the anti-symmetry condition iii). 

In the first model {£,k}k=i 1S a uniformly distributed 
random sequence between [0,1/2], so the size of the 



micro-intervals is random and the function is denned as 

F A {u) = u-$ n + F($ n ) if 116 [£„-i, U , (3) 

where F(£„) is the chaotic map function in Eq. eval- 
uated at £„. 

In the second model, the size of the micro-intervals is 
uniform, i.e. |£„ — £ n _i| = 1/N, but the "height" of the 
function is random. That is, 



F A (u) = u - £ n + K, 



ifti€ [&,-!,£„[ , (4) 



where £ = (£ n _i + £n)/2, and the heigth at the mean 
point, A„, is a random variable distributed uniformly in 
the interval [—6,1 + S], with 6 > 0. Figures and [21 
show examples of these two models. Let us note that the 
random variables and A n are quenched variables, i.e. 
in a given realization of a model they depend on the cell 
Ci (so one should properly write and A« ) but they 
do not depend on time. 

As Fig. H illustrates in model J3J each cell is parti- 
tioned in the same number N of randomly chosen micro- 
intervals of mean size A = 1 /N and in each of them the 
slope of the map (microscopic slope) is one. For A — > 
(equivalently N — > oo) we recover the chaotic system 
but the limit Fa — > F has to be carefully interpreted. 
This kind of modification of the original chaotic system 
is somehow equivalent to replacing circular by polygo- 
nal obstacles in the Lorentz system ^?Jj since the steps 
with unitary slope are the analogous of the flat bound- 
aries of the obstacles. The discontinuities of Fa produce 
a dispersion of trajectories in a way similar to that of a 
vertex of polygons that splits a narrow beam of particles 
hitting on it. However, at variance with the 2-d wind-tree 
model, both randomness and forcing are needed to attain 
a full diffusive behavior. Note that, thanks to the local 
preservation of the antisymmetry with respect to the cell 
cen ter . no net drift is expected and also the Sinai-Golosov 
[271 128L l29l | effect is absent. Also model (@J is non-chaotic 
and the randomness is introduced in the function but not 
in the partition of the unit cells. 

Since Fa has slope 1 almost everywhere, the chaoticity 
condition i) is no more satisfied in these models. How- 
ever, these models can exhibit diffusion provided condi- 
tion i) is replaced by the following two requirements: 

i-a) presence of quenched disorder, 

i-b) presence of a quasi-periodic external perturbation. 

In models (0 and (0J , requirement i-a) is satisfied by the 
randomness in the partition of the micro-intervals, and 
the randomness in the linear function. To fulfill require- 
ment i-b) we introduce a time-dependent perturbation 
into Eq. 



F A (x* 



[x]) + ecos(at) 



(5) 



where e tunes the strength of the external forcing, a de- 
fines its frequency and the index A indicates a specific 
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FIG. 1: Sketch of the random staircase map in the unitary 
cell. The parameter a defining the macroscopic slope is set 
to 0.23. Half domain [0, 1/2] is divided into TV = 12 micro- 
intervals of random size. The map on [1/2, 1] is obtained by 
applying the antisymmetric transformation with respect to 
the center of the cell (1/2, 1/2) 
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FIG. 2: Sketch of model J3| in the unit cell. In this case the 
half domain [0, 1/2] is divided into TV = 12 uniform micro- 
intervals. The map on [1/2, 1] is obtained by applying the 
antisymmetric transformation with respect to the center of 
the cell (1/2, 1/2). In each micro-intervals the map is defined 
by prescribing a linear function with unit slope and random 
height, with 5 — 0.2 (see text). 



quenched disorder realization. We underline again that 
the linear piecewise map Fa explicitly depends on the 
cells Ci due to the disorder, which varies from cell to 
cell. Even though by construction the system is not 
unstable in the Lyapunov sense, it exhibits standard dif- 
fusion. In particular, numerical simulations show a linear 
growth of the mean square displacement. It is important 
to note that, for model (@ a uniform partition of cell 
domains plus the quasi-periodic perturbation is not suf- 
ficient to produce diffusive dynamics. Even in the pres- 
ence of quenched randomness (i.e. random partitions or 
random heights) alone, the system spontaneously settles 
onto periodic or quasi-periodic trajectories and does not 
diffuse. The periodic and quasi-periodic behavior is bro- 
ken down, only when a time dependent perturbation is 



6x10 



4x10 



<(x ) > : 



2x10 



« e = 0.005 
■ £= 0.059 
• e = 0.0034 



0,0 5.0xl0 4 l.OxlO 5 1.5xl0 5 2.0xl0 5 

t 

FIG. 3: Time behavior of ((x*) 2 ) for three different perturba- 
tion amplitudes e, a = 0.23 and N — 100 micro-intervals per 
cell. 



explicitly superimposed onto the quenched randomness. 
Hence, we can say that, in our models, the time plays 
somehow the role of the missin g sp atial dimension with 
respect to the 2-d models of ref . [l7| • We will explore this 
idea in more detail in Sec. V. 



III. DIFFUSION WITHOUT CHAOS, 
NUMERICAL RESULTS. 

In this section, we characterize the diffusion properties 
of model ®, while in the next, we discuss diffusion on 
model (QJ. 

The diffusion coefficient, 



D 



lim — (fx* 



x ) 2 ) 



(6) 



has been numerically computed from the slope of the 
linear asymptotic behavior of the mean square displace- 
ment. The average has been performed over a large num- 
ber of trajectories starting from random initial condi- 
tions uniformly distributed in cell [0, 1] and the results 
are shown in Fig. [3] for three different perturbation am- 
plitudes e. In all performed simultations we present in 
this section we set a — 1; the influence of a on diffusive 
behaviour is discussed in the following section. 

We see that diffusion is absent for e small enough, thus, 
a transition between a diffusive to non-diffusive regime is 
expected upon decreasing e. When the limit in expres- 
sion © exists, coefficient D can be expressed in terms of 
the velo city correlation function as in Green-Kubo for- 
mula nm 



D 



0\2\ 



fc=i 



{v°v k } 



(7) 



where x 4+1 — x % — v % can be formally regarded as a 
discrete-time velocity. The simplest approximation for 
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FIG. 4: Log-Log plot of the dependence of the diffusion co- 
efficient D on the external forcing strength e. Different data 
relative to a number of cell micro-intervals N = 50, 100 and 
150 are plotted vs the natural scaling variable eN to obtain 
a collapse of the curves. Horizontal line represents the theo- 
retical result JS) holding for the chaotic system 11 1 I2t 



D amounts to assuming that velocity correlations decay 
so fast that only the first term of the above series can be 
retained. In the pure chaotic system (I), this condition 
is often verified, and with the further assumption of uni- 



form distribution for the quantity x* 
the approximated expression 

D(a) 



2(1 + a) 



one obtains 



(8) 



that works rather well in the parameter region we ex- 
plored [2q. In the non-chaotic system (J5J this approx- 
imation is expected to hold only in the large driving 
regime limit, where the stair- wise structure of Fa is hid- 
den by the perturbation effects. 

The role of the external forcing on the activation of 
the diffusion process in the non chaotic system has been 
studied by looking at the behavior of the diffusion co- 
efficient D upon changing the perturbation amplitude e 
and the "discretization level" defined by the number of 
the micro-intervals N in each cell. 

Results of Fig. 0| show that D is significantly different 
from zero only for values e > e c , for which the system 
performs normal diffusion. Below this threshold the val- 
ues of D, when non zero, are so small that we cannot 
properly speak of diffusion. When s increases well above 
e c , D exhibits a saturation toward the value predicted 
by formula (jHJ) (horizontal line) for the chaotic system 
defined through The onset of a diffusive regime 

as a threshold-like phenomenon, with respect to the ex- 
ternal perturbation amplitude, is expected. It can be 
explained by noticing that, due to the staircase nature 
of the system, the perturbation has to exceed the typi- 
cal discontinuity of Fa, to activate the local instability 
which is the first step toward the diffusion. Data collaps- 
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FIG. 5: Time behavior of the variance (Eq. |0J for chaotic 
and non chaotic system with a — 0.23. The discretization 
level is N = 100 for the latter. Full lines refer to the non 
chaotic model and from below correspond to initial values 
<7o = 10~ 12 , 10~ 8 , 10 -6 respectively. For a direct comparison, 
full dots represent a 2 (t) in the chaotic case, with a% = 10~ 12 , 
which clearly exhibits a crossover from the exponential growth 
at small times (governed by Lyapunov exponent) to normal 
diffusion at large times. Dotted lines indicate the predicted 
chaotic growth exp(2Af) and diffusion 2D(a)t (see formula|5J. 

ing in Figs 0] confirms this view, because it is achieved by 
plotting D versus the scaling variable eN. This means 
that, for N intervals on each cell, the typical discontinu- 
ity in the staircase-map is 0(1 /N), then the e-threshold 
is e c ~ 1/N. This behavior is robust and does not de- 
pend on the precise shape of the forcing, as we have ver- 
ified by considering other kind of external perturbations. 
However, as will be discussed in the next section, there 
is some connection among the diffusive properties of the 
system, the periodic orbits, the parameter N, and the 
value of the perturbation frequency a. 

The activation of the instability performed by the per- 
turbation can be studied by considering a well localized 
ensemble of initial conditions. We choose uniform initial 
distribution of TV walkers with erg ~ KT 12 ,lCr 8 ,lCr 6 
and monitored the distribution spreading under the dy- 
namics by following the evolution of its variance 

1 M 

^w = ^Z>'-<s*>) a ) (9) 

i— 1 

Note that, in contrast with the average in Eq. © now 
the average is over an ensemble of initial conditions, uni- 
formly distributed on a small size interval. The spreading 
of the initial cloud of points is better understood when 
we make a comparison between the chaotic Q and non 
chaotic system (JjJJ. In the former, the instability makes 
the variance growing exponentially at the rate ln[2(l+a)] 
which is its Lyapunov exponent. This exponential sepa- 
ration lasts until the typical linear behavior of the diffu- 
sion takes place (see Fig. El) - In the non-chaotic system, 



instead, the behavior of cr 2 (t) strongly depends on its 
initial value a 2 and we have to consider two sub-cases, 
whether a > 1/N or a Q <^ 1/N . When ct > 1/N, ba- 
sically one observes the same behavior of the continuous 
limit i.e. a 2 (t) ~ a 2 exp(2Ai) for times t < ~ ln(l/<7o) 
and u 2 (t) ~ t for t > t*. Instead, if <tq <C 1/N we 
can distinguish three different regimes: a) a(t) remains 
constant for a certain time span [0,ti], b) the instabil- 
ity starts being active and the a grows exponentially as 
in the chaotic case, u(t) ~ (l/N)exp[X(t — ti)], c) the 
system eventually reaches the regime of standard diffu- 
sion where a 2 behaves linearly. The crossover time t\ 
between the regimes a) and b) depends on the size of 
the discontinuities of Fa and the specific realization of 
the randomness and, of course, it decreases when either 
N or (To increases. We stress that the above scenario is 
due to the presence of finite size instabilities that in a 
non chaotic Lorcntz gas would correspond to the defo- 
cusing mechanism of a beam by the vertices of polygonal 
obstacles. 

The onset of standard diffusive behavior in model (J5J 
depends on large-time features of the deterministic dy- 
namics. Indeed the presence of periodic trajectories, 
standing or running |.3('|. have a strong influence on the 
diffusion, because the formers tend to suppress diffusion 
while the latter induce a "trivial" ballistic behavior. In 
Ref. [T3, it has been conjectured that the main difference 
between a fully chaotic deterministic system and a non 
chaotic system exhibiting diffusion, as a consequence of 
quenched randomness, is in the different effects the pe- 
riodic trajectories have on the systems. An escape rate 
method 0, El can be applied to analyze these effects. 
This method consists on starting with an ensemble of 
W(0) trajectories initially localized in the region 72. l be- 
tween —L/2 and L/2 and computing the number W(t) 
of them which never left TZl until the time t. For a pure 
random walk process one would expect an exponential 
decay 

W(t) ~ exp(-Ti) (10) 

where the rate 7 is proportional to the diffusion coeffi- 
cient D through the relation 7 cx D/L 2 . In this sense the 
method is an alternative way to estimate the diffusion co- 
efficient from measurement of 7. However possible devia- 
tions from IjlOp generally indicate a non diffusive behavior 
related to the presence of periodic orbits. As clearly seen 
in fig [f)| on all the range of sizes L we have considered, 
we do not observe any substantial deviation from the ex- 
ponential decay that can be associated to periodic orbit 
effects. One can wonder whether the diffusive behavior 
could be merely due to a possible non zero spatial en- 
tropy per unit length of the quenched randomness. To 
unravel the role of static randomness as possible source 
of entropy and thus diffusion, following Ref. [13, we have 
considered a system were the same configuration of the 
disorder is repeated every M cells (i.e. = S,n +M ^), 
so the entropy for unit length is surely zero. This point 
is crucial, because, someone can argue that a determin- 
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FIG. 6: Exponential decay with time of the fraction W(t) /Wo 
of walkers that, at time t, never left the region of length L = 
40, 60, 80, 100, centered around the origin of the lattice of cells. 
Model parameters are a = 0.23, e — 0.1 and N = 100. The 
disorder configuration is replicated each M = 20 cells. The 
walker are initially (t = 0) placed inside the cell [0, 1] 
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FIG. 7: Construction of space-time billiard from a piece-wise 
linear with unit slope map. 



istic infinite system with spatial randomness can be in- 
terpreted as an effective stochastic system, but this is a 
"matter of taste" . Anyway, our system with a spatially 
periodic randonness is deterministic from any point of 
view. 

Looking at the diffusion of an ensemble of walkers we 
obtain that: 

a) there is weak average drift V, that vanishes approxi- 

mately as V ~ 1/M when M goes to infinity. 

b) there is still a diffusive behavior i.e. ((x* — Vt) 2 ) ~ 

2Dt, with a coefficient D very close to the value 
obtained by formula JSJ. 

Let us note that this behavior is obtained by considering 
single realizations of the randomness: the value of the 
drift changes with the realization, in such a way that 
an average on the disorder would give a null drift. The 
fact that V ~ 1/M can be easily interpreted as a self- 
averaging property, for large M. On the other hand if we 
write ((x*) 2 ) = (Vt) 2 + 2Dt, we can define a crossover 
time r c ~ DM 2 , from pure diffusive to a ballistic regime. 
When M is finite but large enough the crossover time 
becomes very large. Therefore, we can safely conclude 
that, even with a vanishing (spatial) entropy density of 
randomness, at least for sufficiently large M, the system 
is able to effectively perform standard diffusion for a very 
long time. However we stress that looking at ((a;* — Vt) 2 ) 
we do not observe any crossover. 

IV. SPACE-TIME BILLIARDS 

As mentioned previously, in the one-dimensional maps 
proposed in this paper, time plays the role of the sec- 



ond dimension in the two-dimensional billiard models 
(e.g., Lorentz and Ehrenfest models). Here we show this 
explicitly by describing how the one-dimensional time- 
dependent maps define billiards in space-time. For con- 
venience, we carry out the construction for model l|4*|l 
where the partition of the unit cells is uniform. However, 
the same construction can be done for model PJl. 

Let [£n— 1> be a micro-interval of a given unit cell. 
Then, for any initial condition xq G [£ n _ i>£n]> the veloc- 
ity of the particle at time t is 

h* n = x t+1 - x l = A„ + e cos(crf) , (11) 

where h n is the height (including the time-dependent per- 
turbation) of the piece-wise linear segment defining the 
map according to Eq. Q. As shown in Fig. [7] in a space- 
time diagram the "world-line" of this orbit is the straight 
line joining the points (x*,t) = (x t+1 , t+1) with slope h l n . 
Accordingly, as shown in Fig. [21 the space-time can be 
decomposed into discrete boxes [£ n -ij£n] x [t,t+ 1], and 
the world-line of successive iterates of the map appears 
as a zig-zag trajectory joining different boxes. The rule 
of the billiard is that a particle entering box n will leave 
this box with speed /i^. That is, in going from box m 
to box n a particle suffers a space-time scattering event 
Sv = h^ 1 —h>m, where Sv = v t+1 — u* is the change in the 
particle velocity. In the space-time diagram the scatter- 
ing rule is represented by drawing at each box a line seg- 
ment with slope given by h* n that according to l|ll(l has a 
random component A„, and a time-periodic component. 

According to this, the time-dependent, zero-Lyapunov 
exponent (slope one) map is equivalent to a space-time 
billiard with scattering rule h^. Note that these mod- 
els are simpler than the usual 2-d billiards (e.g. Lorentz 
gas) because the scattering rule is independent of the an- 
gle of incidence. Although this is just a reformulation of 
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FIG. 8: Space-time billiards without time dependence (e 
0), and with time dependence (e =^ 0). 
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FIG. 9: Periodic orbits in a space-time billiard with rational 
perturbation frequency 



the same model, the billiard construction provides use- 
ful insights into the problem and allows the connection 
with the two-dimensional billiard problems. When the 
map is time-independent (e = 0) the scattering rule de- 
pends only on n {h^ — h q n for any integers p and q) and 
as shown in Fig. [HJ- (a) the billiard structure is simply a 
copy of the billiard at t = 0. In this case, despite the 



spatial randomness, the lack of a time dependence leads 
to a highly structured space-time billiard in which peri- 
odic orbits are ubiquitous. As a consequence, in this case 
there is no diffusion, a result already mentioned before. 
On the other hand, as Figure IHKb) shows, when e ^ 
the billiard exhibits space-time disorder which leads to 
the possibility of diffusive behavior. 

As discussed in Ref. 01 f° r polygonal two-dimensional 
billiards, periodic orbits play an important role. The 
existence of periodic orbits in the space-time billiard is 
intimately connected with the frequency a of the time 
periodic perturbation. If {x*} is a periodic orbit of period 
M [x t+M — x l for any t), the billiard must be g-periodic, 
that is ft-^ +<? = h l n for any t, where M = mq and m is a 
positive integer. According to Eq. this will be the 
case provided 



a = 2 it v . 



(12) 



with v = p/q. Figurc[5]show two examples of periodic or- 
bits with e = 0.2. In panel (a), a — 2tt(3/A), and the pe- 
riod of the billiard is the same as the period of the orbit, 
q = M = 4. In panel (b), a = 27r(2/5), and the period 
of the orbit is twice period of the billiard M = 2q = 10. 
Thus, for rational v the billiard shows a periodic struc- 
ture in time, and periodic orbits are likely to exist leading 
to the suppression of diffusion. On the other hand, when 
v is an irrational number the billiard structure never re- 
peats exactly in time and strictly speaking there are no 
periodic orbits. Motivated by this, it is possible to re- 
lated the degree of temporal disorder of the billiard to 
the degree of irrationality of v as determined, for exam- 
ple, by the continued fraction approximation of v. In 
this regard, billiards with irrational frequencies v which 
are easily approximated by rationals are "less disordered" 
than billiards with hard to approximate frequencies. To 
explore this ideas Fig. EH shows the time evolution of 
the variance for different values of the perturbation fre- 
quency in the same random realization of model I0J with 
N — 40 and e = 0.4. For the perturbation frequencies, 
a = (p/q)2ir, we took the continued fraction expansion 
of the inverse golden mean 



- v= 1/2, 2/3, 3/5, ...F n ^/F n 



V5-l)/2 = 7 - 1 
(13) 

where Fo = F\ = 1, and F n = F n -i + F„-2- The chose 
of the golden mean is motivated by the fact that this 
number is the hardest to approximate with rational. In 
the model the disorder configuration is repeated every 
M cells, and for small M the recurrence time of an or- 
bit (i.e. the time to return to the same microinterval, 
module M) decreases. Since we are interested in the 
role of periodic orbits, here we consider the limit case 
M = 1 . That is, in the simulations reported in Fig. 1101 
the disorder configuration in the unit interval [0, 1] was 
copied to all the cells [£,£ +1]. As expected the figure 
shows the suppression of diffusion for rational frequen- 
cies. However, more interesting is the fact that there is 
a tendency for the onset of the suppression of diffusion 
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FIG. 10: Mean square displacement as function of time in 
model (2J for a sequence of values of the perturbation fre- 
quency q approaching the inverse golden mean I/7. The 
heights of map are uniformely distributed in [0, 1] (5 = 
see text). 

to increase as v approaches 1/7. This is related to the 
fact that for rationals with lower denominators the time 
period of the billiard is small, and the probability that 
the recurrence time of an orbit coincides with a multiple 
of the period of the billiard is high. The recurrence time 
is also related to the value of TV. When TV is large, the 
probability for an initial condition to return to the same 
micro-interval (module M) is small, and periodic orbits 
are more scarce. This explains why for large enough TV 
and M, diffusion might not be suppressed even for ratio- 
nal frequencies with large denominators. This suggests 
that, for large enough TV and M, deviations from diffu- 
sion are unlikely to be observed numerically in the case 
of rational frequencies with large denominators. 

V. CONCLUSIONS AND REMARKS 

From the above results we can conclude that two basic 
ingredients are relevant for diffusion, 1) a finite size in- 
stability mechanism: in the chaotic system this is given 



by the Lyapunov instability, in the "stepwise" system 
this effect stems from the jumps 2) a mechanism to sup- 
press periodic orbits and therefore to allow for a diffu- 
sion at large scales. In the presence of "strong chaos" 
(i.e. when all the periodic orbits are unstable) point 2) 
is automatically guaranteed thus the periodic orbits for- 
malism |y, 0, LJ] can be fully applied to compute the dif- 
fusion coefficient. In non chaotic systems, however, like 
the present one, the mechanism 2) is the result of a com- 
bined effect of quenched randomness and quasi-periodic 
forcing. 

In summary, we have introduced and studied mod- 
els exhibiting diffusion in the absence of any source of 
chaotic behavior. The models represent the 1-d analogue 
of non chaotic Lorentz gas (i.e. with polygonal obstacles) 
discussed by other authors in connection with the debate 
around experimental evidences for a distinction between 
chaotic and stochastic diffusion. Our results indicate that 
when the chaotic instability condition (positive Lyapunov 
exponent) is replaced by the presence of finite size insta- 
bility (non positive Lyapunov exponent) we need both 
an external quasi-periodic perturbation and disorder for 
preventing the system falling into either a periodic/quasi- 
periodic evolution or a drift dominated behavior. 

We believe that this kind of behaviour is rather gen- 
eral, in the sense that chaos does not seem to be a neces- 
sary condition for the validity of some statistical features. 
This is in agreement with the recent results of Ref. [2^| , 
where the applicability of the generalized Gallavotti- 
Cohen fluctuation formula |3lT ] has been proven for non 
chaotic systems too. 
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